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Abstract: VERITAS is an array of four identical télescopes designed for detecting and 
measuring astrophysical gamma rays with énergies in excess of 100 GeV. Each télescope uses 
a 12 m diameter reflector to collect Cherenkov light from air showers initiated by incident 
gamma rays and direct it onto a 'caméra' comprising 499 photomultiplier tubes read out 
by flash ADCs. We describe here calibration methods used for determining the values of 
the parameters which are necessary for converting the digitized PMT puises to gamma-ray 
énergies and directions. Use of laser puises to détermine and monitor PMT gains is discussed, 
as are measurements of the absolute throughput of the télescopes using muon rings. 



Introduction 

Like ail gamma-ray detectors which use the 
atmospheric Cherenkov technique, the VER- 
ITAS instrument is fundamentally quite sim- 
ple. Each of its four télescopes consists of a 12 
m reflector which directs Cherenkov light from 
air showers onto a matrix of 499 photomulti- 
plier tubes (PMTs) which are read out using 
500 MSample/s flash-analog-to-digital convert- 
ers (FADCs). In order to translate the digi- 
tal information emerging from the FADCs into 
a form which can be used to select gamma- 
initiated showers from background and déter- 
mine the energy and direction of the inci- 
dent gamma ray, one needs calibration con- 
stants. Thèse 'constants' (which are not, 
strictly speaking, constant) need to be deter- 
mined when commissioning the detector and 
monitored and adjusted periodically during the 
lifetime of the project. In this paper we de- 
scribe techniques employed by the VERITAS 
collaboration to accomplish this task; two tech- 
niques use a laser to détermine the absolute 
gains of the PMTs and one uses Cherenkov 
images, gêner ated by isolated muons, for inter- 
telescope calibration and détermination of ab- 
solute throughput. Some of thèse issues are ad- 



dressed independently with a remote LIDAR- 
like System described elsewhere in thèse pro- 
ceedings [1]. 

The VERITAS Laser System 

For flat-fielding and gain monitoring, VERI- 
TAS uses a nitrogen laser (A = 337 nm, puise 
energy 300 /iJ, puise length 4 ns). The beam 
is sent through neutral density filters arranged 
in two sequential wheels, with 6 filters each, 
such that transmissions ranging from less than 
0.02% to 100% may be chosen. It is then di- 
vided, approximately equally, among 10 optical 
fibres, four of which are routed to opal diffusers 
located on the optical axes of the télescopes, 4 
mètres from the PMTs in the caméras. A fifth 
fibre supplies light to a PIN photodiode to pro- 
vide a fast external trigger for FADC readout; 
self-triggers using only PMT information are 
also used for some applications. There is, at 
présent, no independent monitor for measur- 
ing the pulse-to-pulse fluctuations in the laser 
intensity (typically 10%); thèse are monitored 
using a sum over a large number of PMTs in 
each caméra. 
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Nightly Laser Runs 

A five-minute, 10 Hz laser run at nominal in- 
tensity is taken at the beginning of each observ- 
ing night. The data obtained are used primar- 
ily for monitoring gain évolution and checking 
for problems. Other tests, described below, 
are done less frequently. Since the opal dif- 
fuser spreads the laser light uniformly (to bet- 
ter than 1%) over the face of the caméra, the 
puises can be used for flat-fielding the response 
of the channels. The high voltages of the indi- 
vidual PMTs are adjusted so that the average 
puise size in each channel is the same for ail 
channels. A PMT's average puise size dépends 
on the product of its photocathode's quantum 
efflciency and the efflciency for photoelectrons 
to be collected by its first dynode, as well as 
on the gain in the électron multiplier stage. To 
a lesser extent it dépends on the reflectivity of 
the Winston-cone light concentrator in front 
of each PMT. The average puise sizes are cal- 
culated and written to a database for use in 
off-line analysis. 

The gain of the électron multiplier can be 
tracked separately using the daily laser data 
using the method of photostatistics. In this 
method, we remove laser fluctuations using 
a sum-over-PMTs monitor and the effects of 
electronics noise and night sky background are 
measured in runs with zéro laser intensity and 
unfolded. Then, to first order, we can state 
that the mean charge in a laser puise is given 
by /i = GN pe with G the gain and N pe the 
mean number of photoelectrons arriving at the 
first dynode. Assuming that only Poisson fluc- 
tuations in N pe détermine the width, 
the charge distribution, we have a - 
Thus we can solve for gain as G = a 1 j \i. Tak- 
ing into account statistics at the other dynodes, 
which are in gênerai described by a Polya dis- 
tribution, leads to a correction factor such that 
G = a 2 /fi/(l + ot 2 ) where a is the width pa- 
rameter which would resuit from injecting only 
single photoelectrons into the dynode chain. 
For our PMTs and their associated dynode 
voltages we simulate a = 0.47, which results 
in a revised estimate for multiplier gain of 
G' = 0.82 cr 2 //i. 



As a check on this model, note that /i = G' N pe 
or N pe = fi/G' . This quantity is plotted for a 
représentative PMT, in figure 1, as a function 
of applied high voltages in steps from nominal 
HV. Except perhaps for an effect due to in- 
creased first dynode collection efflciency due to 
increased HV, we do not expect N pe to change 
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Figure 1: Mean number of photoelectrons cap- 
tured by the first dynode of a représentative 
PMT vs déviation from its nominal high volt- 
age setting. The flat line is to guide the eye and 
emphasize that there is no significant change. 



Single Photoelectrons 

An alternative method for determining PMT 
gain is to directly measure the position of the 
single photoelectron peak in a puise size spec- 
trum. Again, this gives the gain of the électron 
multiplier structure (and any downstream elec- 
tronics) and does not include effects of the pho- 
tocathode. To résolve the single photoelectron 
peak, we take spécial laser runs at very low 
intensity where the average number of photo- 
electrons resulting from each laser puise is less 
than 1.0. The resulting spectrum consists of 
a pedestal, the single photoelectron peak, and 
small admixtures of two, three etc peaks with 
the relative sizes of each component prescribed 
by Poisson statistics. We also rely on the 
constraint that the multi-photoelectron signais 
can be fit with the same parameters (mean and 
width) as the single photoelectron peak (up to 
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multiplicative factors). Allowing only a small 
number of free parameters results in robust fit s 
to the spectra, an example of which is shown 
in figure 2. In this example the relative width 
of the single photoelectron peak is found to 
be 0.48. The average from a group of chan- 
nels is 0.47. The data for this study were ob- 
tained using twice the normal gain, in order to 
résolve more clearly the single photoelectron 
peaks. With such a gain the simulation pre- 
dicts a relative width of 0.44. 




Integrated Charge [d.c] 

Figure 2: A puise size spectrum made with 
highly attenuated laser puises and raised high- 
voltage. The single photoelectron peak is 
clearly visible as the structure next to the 
pedestal, which is the dominant feature. The 
data are fit with a sum of Gaussians as de- 
scribed in the text. 

In order to maintain good signal-to-noise for 
this measurement, we cover the caméra with a 
thin aluminum plate with a 3 mm hole drilled 
at the location of the centre of each PMT. This 
reduces the night sky background to the point 
where it is negligable compared with the laser 
light. Indeed, with the télescope in stow po- 
sition, one can perform single photoelectron 
laser runs in the présence of moonlight. This 
is an important considération given that suf- 
ficient statistics (^50000 shots) require nearly 
an hour of running. 

The gain values which resuit from the method 
of photostatistics and from the single photo- 
electron fitting are in units of digital counts 
per photoelectron. A comparison of the two 




Laser 



Figure 3: A comparison of gains determined 
using photostatistics (abscissa) with those de- 
termined from single photoelectron fitting (or- 
dinate). The slope of the corrélation is approx- 
imately 1.1. It should be 1.0; the discrepancy is 
an indication of the présent scale of the system- 
atic error of the gain-measuring procédures. 

methods is shown in figure 3 where results from 
télescope 1 are shown. 

The data points in this figure highlight the dif- 
férence between the 'multipler gain', which in- 
cludes everything starting from the first dyn- 
ode, and the 'overall gain' which also includes 
the light concentrator cônes and the photo- 
cathodes. Since the PMTs have ail been flat- 
fielded according to the overall gain, the dis- 
persion seen along the corrélation line in fig- 
ure 3 is due to channel-to-channel différences 
in thèse 'front end' components. 

Muon Rings 

Local muons are normally a nuisance for 
Cherenkov télescopes but they can be use- 
ful in providing a measurement of the opti- 
cal throughput of the detector [3, 2]. Muons 
passing through the centre of the télescope 
with trajectories parallel to its optical axis will 
produce azimuthally symmetric rings in the 
caméra. The rings will have radii given by 
the Cherenkov angle of the muons (maximum 
value about 1.3 degrees) and the total number 
of photons expected in the ring can be calcu- 
lated from the measured value of this angle. 
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Figure 4: A muon image recorded by VERI- 
TAS télescope 1. Puise size in each PMT is 
colour-coded - the azimuthal nonuniformity is 
the resuit of a non-zero impact parameter. 

Muons with non-zero impact parameters will 
produce arcs with an azimuthally dépendent 
photon density and muons arriving at an angle 
with respect to the telescope's axis will give 
rise to arcs with centres that are offset from 
the centre of the caméra. 

Muon ring images can be obtained from nor- 
mal data where they occur as part of hadronic 
showers. The images are cleaned (channels are 
required to have a minimum puise size and to 
be next to other channels with non-zero charge, 
otherwise they are set to zéro) and a ring is fit 
to the image. Further cleaning of the images, 
where charge deposits far from the fitted ring 
are suppressed, removes light from other com- 
ponents of the shower of which the muon was 
a member. After this second cleaning the ring 
parameters are re-calculated. A cleaned image 
of a complète muon ring is shown in figure 4. 

Since the morphology and location of the muon 
ring allow the muon's trajectory to be calcu- 
lated, it is possible to predict with précision 
the number of Cherenkov photons that should 
be collected by the caméra. This requires 
knowing the reflectivity of the mirror facets, 
shadowing effects due to the caméra support 
structure, etc so the detector response to local 
muons is a good check on our understanding of 
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Figure 5: Detected charge in muon arcs, nor- 
malized to their lengths, for VERITAS téle- 
scopes 1 (histogram) and 2 (data points), 
showing that they are well matched. 

the instrument. Absolute calculations are still 
in progress but certain relative measurements 
have already been implemented, such as inter- 
telescope calibration and month-to-month sta- 
bility checks. An example is shown in figure 5 
where we histogram the summed charge in each 
muon arc, divided by its length, for two téle- 
scopes in the array. The overlap of the two his- 
tograms, normalized by the number of entries, 
shows that the télescopes are well balanced. 
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